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matrices and log-gas systems at infinite density 
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Abstract 

It is shown how the universal correlation function of Brezin and Zee, and Beenakker, for 
random matrix ensembles of Wigner-Dyson type with density support on a finite interval can 
be derived using a linear response argument and macroscopic electrostatics. The analogous 
formula for ensembles of unitary random matrices is derived using the same method, and the 

o ' 

fi , corresponding global fluctuation formula for the variance of a linear statistic is presented. 



The result for the universal correlation is checked using an exact result for the finite system 
two-point correlation in the Dyson circular ensemble at all even j3. 

PACS numbers: 72.10.Bg, 05.40.+J, 05.60,+w 

1 Introduction 

A large portion of random matrix theory is concerned with the eigenvalue probability density 

function 

N 

Y[ e -W(*i) J] \xk-xjf, XjCR (1.1) 

1=1 l<j<k<N 

where (3 = 1,2 or 4 according to the matrices being real symmetric, Hermitian or quaternion 
real respectively (see e.g. [1]). An alternative interpretation of (1.1) is as the Boltzmann factor 



* Supported by the Australian Research Council 



of a one-component log-potential Coulomb system confined to a line with potential energy 



N 



J2 v ( x i)- Yl l °s\ x k 



X2) 



1=1 l<j<k<N 

The two body term in (1.2) is the electrostatic energy of the Coulomb repulsion between like 
charges, while the one body potential acts to confine the particles to some interval of the real 
line. 

In particular, consider the situation in which V{x) is infinite outside an interval (a, b) so that 



Xj G (a, b) 



(1.3) 



This setting has recently been discussed by Brezin and Zee [2], in the special case (3 = 2 and 
V{x) an even polynomial about the centre of the interval (see also [3]), and by Beenakker [4], 
for general (3 and V(x). These authors have found that the two-point correlation, p T (x,x') say, 
satisfies the universal formula 



lim 

N—*oo b — a 



N rx'+(b-a)/N 



dtpjf(x,t) 
1 (a + b)(x + x')/2 — ab — xx' 



x / x' 



(1.4) 



(3tt 2 (x - x') 2 [(x - a){b - x)(x' - a){b - x')] 1 / 2 ' 
Remark: The function p T (x,x') above equals —T2(x,x') in the notation of Beenakker [4]. 

The l.h.s. of (1.4) corresponds to the two point correlation in the infinite density limit, smoothed 

over the average interparticle spacing. 

As a consequence of (1.4), Beenakker [4] has derived the global fluctuation formula 

N 



lim var 



£/(; 



i=i 



(3lT 



1 f b f b i 
VI dx I dx 



{x' — a)(b — x') 



l2/W ^/M 



x — x' dx 1 



(1.5) 



(x — a)(b — x) 
where V denotes the principal part, for the variance of a linear statistic. 

The study of random matrix ensembles with density support confined to a finite interval 
follows the study [5-9] of ensembles with density support on the half interval [0, oo), which are 
of relevance to fluctuation phenomena in mesoscopic systems. In this setting it is not necessary 
to consider an infinite density limit. Rather, a scale can be chosen so that the number of 
eigenvalues near the "edge" x = remains finite [7] in the N — » oo limit. The asymptotic 
expansion of the corresponding two-point correlation satisfies a formula analogous to (1.3): 



p T (x,x') 



1 



x + x' 



2TT 2 f3Vxx 7 (x-x') 2 



(1.6) 



for the leading non-oscillatory behaviour. Beenakker [9] (see also [10]) has used (1.5) to derive 
the analogue of (1.4) for the fluctuation of a linear statistic. 

Recently [11] a new derivation of (1.5) has been given, applicable to all log-potential Coulomb 
systems in their conductive phase, which is based on macroscopic electrostatics and linear re- 
sponse theory. It is our purpose herein to first rederive (1.4) using the method of [11], and then 
to consider the analogues of (1.4) and (1.5) for unitary random matrices (or equivalently log-gas 
systems confined to a circle). For the one-component log-gas on a circle at even /3, the two point 
correlation has been evaluated exactly [12] in terms of a /^-dimensional integral for all particle 
numbers N, which allows an explicit verification of the analogue of (1.4). 

2 Macroscopic electrostatic argument 

Our objective in this section is to derive (1.4). For simplicity of notation we consider the 
symmetric case a = —b and scale the variables x i— ► bx and x' h^ bx' , so that the spectrum has 
support in (—1,1). The r.h.s. of (1.4) then reads 

1 1 — xx' / 

/37r 2 (x-x / ) 2 [(l- x 2)( 1 _ x /2)]|' 

There is no loss of generality is specializing to this case as the r.h.s. of (1.4) can be reclaimed 
from (2.1) by rescaling the variables x t— ► 2{x — (a + b)/2)/(b — a) (and similarly x'). Rather 
than restricting ourselves to systems with probability density function (1.1) we will consider any 
log-potential Coulomb system confined to the interval (—1,1) which consists of point particles 
and is in a conductive phase. 

Our starting point is the linear response argument of [11], which says that if the system is 
perturbed by adding an external charge 5q at a point r', then 

< $(r) > - < $(r) > = -P6q < $(f)$(P) > T (2.2) 

where <3?(r) is the potential at the point r due to the induced charge distribution, <>o refers to 
the value of the average before the perturbation and < AB > T =< AB >o — < A >o< B >$. 
To obtain from (2.2) the charge density- charge density correlation 

< <j{x)o{x') > T (2.3) 

we hypothesize that in the high density TV — ► oo limit, the system behaves as a conducting 



interval from —1 to 1 along the x-axis, obeying the laws of two-dimensional macroscopic elec- 
trostatics. Now, for a point x on the x-axis, macroscopic electrostatics gives 

E+(x)-E y (x) = 2Tra(x), (2.4) 

where the y-components of the electric fields are given in terms of the potential by 







d 



E vW ■■= 7r-*(0 „_ n+ E-(z) := j-m __ (2.5) 



dy y=o+ y dy 



y=0- 



Hence 



-/3Sq < a(x)a(x') > T 
= -P^J^L (< E^{x)E+{x>) > T -< E+(x)E-,(x') > T 
2 ( d d < *,-> */-. 



(2vr) 2 \dydy' 

!:^7(<^)>-< <1> ('0>o) „ n ) (2.0) 



9y 5y 
where to obtain the first line we have used the symmetry 



y=0+,y'=0+ 
y =0+,y'=0- 



< E+(x)E*,(x>) > T =< E-{x)E-°(x') > T 

for s = +, — and to obtain the second line (2.5) and (2.2) have been used. 
By the hypothesis of the applicability of macroscopic electrostatics, 

< $(f) > - < $(f) > 

is simply the potential at a point r due to the charge induced on the conducting interval 
(—1,1) by the external charge 5q. The calculation of this potential is a standard problem in 2d 
electrostatics [13, Chapter 4, exercise 5]. One first considers the potential at r, \f\ > 1, due to a 
unit charge at r', \r'\ > 1, and a conducting unit disk centred at the origin. By the method of 
images, this potential is 

$(z,z') = -log(\z' - z\/\z' - l/z*\) - log \z'\ (2.7) 

where z{z') is the complex coordinate corresponding to r = (x,y)(r' = (x',y')), and it is 
assumed the conductor carries no net charge. Then one maps the exterior of the unit disk to 
the plane, with a cut from —1 to 1 along the x-axis by the conformal mapping 

w = z +- and thus z = w + (w 2 - 1) 1/2 , (2.8) 



where the principal branch of the square root is chosen. 

Using (2.8) to substitute in (2.7) for z (and similarly z') we obtain the potential at r due to 
the charge induced on the conducting interval and due to the external charge at r' . Subtracting 
the potential due to the external charge we therefore have 

< $(r) > - < $(f) > 
= Sq(\og\w - w'\ -log|u/-w + (u/ 2 -l) 1/2 -(u; 2 -l) 1/2 | 

+ log \w' -w* + (w' 2 - 1) 1/2 + ((w 2 - 1) 1/2 )*| - log \w' + (w' 2 - 1) 1/2 |) 



where 



w = x + iy and vJ = x + iy' 



and we have used the formula 



l/( U ; + ( U ; 2 -l)V2 )=u; _ (u; 2_ 1) l/2 



(2.9) 



(2.10) 



(2.11) 



Substituting (2.9) in (2.6), it remains to perform the derivatives. Now, for —1 < x < 1 and 
0<y< 1 



(w 2 - I) 1 ' 2 

so differentiating (2.9) with respect to y gives 




x z -l 1/z ie 



l/2-ixy/(x 2 -l) 



(2.12) 



dy 



( < $(r) > - < 3>(r) > 



J / =0^ 



Sq -Re— 



w — x 



|l_ x 2|l/2 



Re 



u 



W' + (W> 2 - 1)V2 - u 



where 



u := x + i\x — l\ ' 



(2.13) 



(2.14) 



Using (2.11) with w,x,y replaced by w',x',y', (2.12) then gives 



~{ < *(f) > - < m >o 



Sq 



1 



+ 



y=0+,y'=0+ 

2 



( X '-X)2 |1_ x 2|l/2|! _ x /2|l/2 



Re 



(u' - n) z 



where 



n'^x' + ilx' 2 -!! 1 / 2 
Also, since for — 1 < x' < 1 and — 1 <C y < 

{w' 2 - l) 1 / 2 



^ -iZ-lli/VAV^-i), 



(2.15) 



(2.16) 



(2.17) 



we have 



d_d_ 

dy dy 



-{ < $(f) > - < $(r) > ) 



Sq 



1 



y =0+,y'=0- 

2 



-Re 



/* 
u u 



(u'* — -u) 5 



(2.18) 



(x'-x) 2 |1 -x 2 \V 2 \l- x >2 \V 2 ' 
We now substitute the evaluated derivatives (2.15) and (2.18) in (2.6) to obtain for the 
charge-charge correlation 



< a(x)a(x') > J 



1 



r./2 1 1/2 



Re 



u'u 



+ 



/* 
u u 



(u' — u) 2 (u 1 * — u) 



(2.19) 



pir 2 \l - x 2 \V 2 \l 

Straightforward simplification of the r.h.s. of (2.19) using the definitions (2.14) and (2.16) of u 
and v! reclaims (2.1). 



3 Ensembles of unitary random matrices 

The Dyson circular ensembles of unitary random matrices (see e.g. [1]) have eigenvalue proba- 
bility density function 

TJ | e ^_ e i0 3 |/3 (3.1) 

l<j<k<N 

where (3 = 1,2 or 4 according to the ensemble being invariant under orthogonal, unitary or 
symplectic similarity transformations respectively. The probability density is identical, up to 
a multiplicative constant, to the Boltzmann factor of a one-component log-potential Coulomb 
system confined to a unit circle. It is particularly simple to apply the argument of Section 2 to 
this setting and thus deduce that the corresponding two-point correlation p T (9,9') satisfies the 
universal formula 



N t-e'+^/N 
hm — / d<pp N {9,(j)) 

iV^oo 2ir Jo 



(3.2) 



/3{2tt) 2 sin 2 ((0 - 6')/2) 
Let us briefly indicate the required working. 

First, since the domain is a circle, in polar coordinates the perpendicular component of the 
electric field is the r-component, and thus the analogue of (2.4) is 



E+(e)-E-(9) = 2va(d) 



(3.3) 



Let us now assume that the source point r' is outside the unit circle. Then, if the field point r 
is also outside the circle, the total potential created is the same as that with the circle replaced 



by a unit conducting disk and is thus given by (2.7). Subtracting the charge-charge portion of 
this potential then gives 

< $(r) > - < $(r) > =<5g(log|z'-l/z*| -log|z'|) (3.4) 

For the field point inside the unit circle, there is no net potential due to the source point outside 
the circle. Subtracting the charge-charge potential gives in this case 



< <&(r) > - < $(f) > = 5qlog\z' - z\ 



(3.5) 



The analogue of (2.6) can now easily be evaluated to derive the r.h.s. of (3.2). 

The result (3.2) can be used to derive for random unitary matrices (or equivalently, log- 
potential Coulomb systems on a circle) the analogue of the global fluctuation formula (2.6) . 
For general systems of point particles on a unit circle 

v 



lim var 

JV— >oo 



EM 

3=1 

2tj- p2n 

liin r dd d6'f(e)f(e')(p%(e-6')+ P 5(sm{6-e')/2) 



.. (3.6) 

N^odJq Jo \ I 

where p := N/2n and the delta denotes the Dirac delta function. Introducing the Fourier series 

oo 

i{e-6')n 



lim (p 1 N (e-6') + p6(8m(6-e')/2) 



E 



Sn.C 



this reads 



where 



lim var 

N— >oo 



N 

EM 



/ j s n\Jn 



(3.7) 



(3.8) 



And 



(3.9) 



/„:= r def(ey 

Jo 
The coefficients s n are given by the usual inversion formula of (3.7). In the inversion formula 

we replace p (9) by its smoothed value (3.2) written as a second derivative: 

1 d 2 



p'iO) 



log|sin0/2|, 0^0 



ir 2 [3 d0 2 

We then calculate the coefficients by a formal differentiation by parts: 

1 1 



(3.10) 



tt 2 /3 2W-^ W 



/ d 2 
d0( — log|sin0/2|)e 



iriO 



1 



7T 2 /T 

\n\ 
2vr 2 /3 



-»-)— / d6 (log|sin0/2|)e 



2n J-k 



inO 



(3.11) 



Hence, after substituting (3.11) in (3.8), we obtain the desired global fluctuation formula 

1 



lim var 



N 

Em) 



^2g E Nl/nl 2 (3-12) 



A recent exact calculation [14] allows the formula (3.2) to be illustrated for the one-component 
log-gas on a unit circle at all even values of (3. The exact calculation gives that the two particle 
distribution function Pjv(0, 6) for this system is equal to a /^-dimensional integral [14, eq. (3.15)|^| 
with 7 = /?]: 

Pn(0,9) = N{N L ; l) (2sin^/L)^/3)e-^ (JV - 2)/L 



r p 

x y^ i ^dh . . . dtp J] [1 - (1 - e 2 ^/ L )t i ] 7V - 2 D 2//3)2//3>4//3 (t 1 , ...,t p ) (3.13) 

where 

D 2 /i3,2/PA/0(ti,...,tp) := UKl-tjfa] 1 - 2 '? II l<*-*il 4/ ^ ( 3 - 14 ) 

3=1 l<j<fe</3 

and L = 27r. The exact value of the normalization Bjsr(/3), which is given explicitly in [14], is 
not required for our discussion below. 

We seek to expand (3.13) for large- N so that the integral on the l.h.s of (3.2) can be evaluated. 
In the large- A r limit we see that the iV-dependent factors of the integrand in (3.13) take on their 
maximum values at tj = and tj = 1. The large- N expansion can be obtained by expanding 
the integrand about the maximum values. Note in particular that for tj ~ and L = 2n 

[1 - (1 - e 2m0 / L )t,] N - 2 ~ e -(.N-W-e»ytj (3 15) 

while for tj ~ 1 and L = 2ir 

[1 - (1 - e 2 ^ L )tj] N - 2 ~ e i(N-2)e e -(N-2)(i-e^)(i-t 3 ) (3 16) 

At this stage we observe that the problem is very similar to one we have solved recently 
[12]. The solved problem is the task of calculating the large- a; asymptotic expansion of the 
thermodynamic limit form of the two-particle distribution. Now in the thermodynamic limit 

N, L — > oo (N/L = rj fixed) (3.13) becomes 

P 



2mtjrjx 



p(0, x) = V 2 (rjxfA(l3)e- m ^ x J fj dh . . .dtp f[ e 2 

x D 2//3,2/l3,4/l3(h, ■■■ ,tp) (3-17) 



lr rhis equation is missing a factor of £ -™vi(n-2)/n . Q f gq ^ ^ rf j 14 j 



from which the large-x asymptotic expansion can be shown to be of the form 



p T (0,x)~r) 2 



1 a^ A > >-^ cos2TTT]xn I / \ bn 



»- ' y^ wa^/nixu (0) u n 

^ lV.^4 "•" ■ • • "*" Z^ („„\4n*IB n ^ („„\2 + 



/3(7r?7x) 2 (r/x) 4 ^ (r/x) 4 ™ 2 //? ^ n (ryx) 



^ sin27r??xn / ( o) , c " } , 



(3.18) 



The method used in [12] to compute (3.18) was to expand the integrand in (3.17) about the 
endpoints tj = 0, 1. We thus see that this problem is almost identical to the task of computing 
the large- TV expansion of (3.13). 

Indeed, by comparing (3.13) and (3.17) and examining (3.15) and (3.16) we see that the 
large- N expansion of (3.13) is identical to the large-x expansion of (3.17) and is thus given by 
(3.18) provided we make the replacements 

vrr/x^ 2iVsin6V2, rj i-> N/2tt. (3.19) 

Substituting the resulting expansion in the l.h.s. of (3.2) we see that the oscillatory terms all 
integrate to zero, and all but the first of the non-oscillatory terms vanish in the N — > oo limit. 
The first non-oscillatory term gives the r.h.s. of (3.2) as required. 
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